Exact infinite-time statistics of the Loschmidt echo for a quantum quench 
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The equilibration dynamics of a closed quantum system is encoded in the long-time distribution 
function of generic observables. In this paper we consider the Loschmidt echo generalized to finite 
temperature, and show that we can obtain an exact expression for its long-time distribution for a 
closed system described by a quantum XY chain following a sudden quench. In the thermodynamic 
limit the logarithm of the Loschmidt echo becomes normally distributed, whereas for small quenches 
in the opposite, quasi-critical regime, the distribution function acquires a universal double-peaked 
form indicating poor equilibration. These findings, obtained by a central limit theorem-type result, 
extend to completely general models in the small-quench regime. 
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Introduction Imagine an isolated quantum system, 
say the laboratory, prepared in a state po. Accord- 
ing to the laws of quantum mechanics, the state will 
evolve unitarily into p (t) . The average result of a mea- 
surement of an observable O will be the time average 
(O (t)) := J^{0 it))dt, where T is the measurement 
time. Since T is much larger than the microscopic time 
scales of the system it is often set to infinity for math- 
ematical clarity. Now, the postulates of statistical me- 
chanics assert that the time-averaged expectation value 
is indistinguishable from that obtained using the statis- 
tical microcanonical ensemble. Although this postulate 
is confirmed by a number of numerical simulations (see 
G-g- [J), to date no explanation exists for why this is 
so. In other words, the mechanisms of thermalization in 
quantum systems are unknown (though there exist pos- 
sible approaches such as the eigenstate thermalization 
hypothesis [2| or normal typicality ^]). 

In such a context it is important to have exact re- 
sults, at least for some particular cases, which can serve 
to guide our intuition. Ideally one is interested in the 
full, long-time statistics of a generic observable {0{t)). 
This article provides a result in this direction. Namely, 
concentrating on the Loschmidt echo, we will obtain its 
exact, long-time distribution function and investigate the 
effects that proximity to critical points has on the equi- 
libration dynamics. In the thermodynamic limit, also 
called the off-critical regime, i.e. when the system size is 
much larger than all length scales of the system, we will 
see that a central limit theorem result applies leading to 
universal Gaussian equilibration. In the opposite regime 
of quasi-criticality, where the correlation length is equal 
to or larger than the system size, we will again find uni- 
versal behavior, although one in which fluctuations are 
large and thermalization does not occur. 

The scenario we consider here is that of a quantum 
quench, generalized to the mixed case. A closed system 
is initialized in the state po commuting with the Hamil- 
tonian Hq. The system is then instantaneously quenched 



and left to evolve according to Hamiltonian Hi. This is 
an important generalization, since in principle there is 
no reason why the "initial" state of the system should be 
pure. In particular, for its experimental relevance we will 
use Gibbs initial states po ~ e~^^°. Such a situation is 
in fact often realized in the laboratory by first thermaliz- 
ing the system by putting it in contact with an external 
reservoir and then detaching the reservoir. 

The quantity we consider is the Loschmidt echo (LE), 
which generalized to the mixed case is given by 

C (t) = F{p (t) , po) , F{p, a) - (tr VT^W^) ' ■ 

Here F is the Uhlmann fidelity Q which characterizes the 
degree of distinguishability between two mixed states. 
Note that if either (or both) of p and a is pure, the 
Uhlmann fidelity simplifies to F{p,a) = tr (per). 

The quantum XY chain The model we investigate 
here is the quantum XY chain in a transverse magnetic 
field, 
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A Jordan- Wigner transformation brings Eq. ([T]) to a 
quadratic form in Fermi operators Ci, and hence can be 
exactly diagonalized. At zero temperature the model ([T]) 
displays two kinds of quantum phase transition lines in 
the {h, 7) plane. For /i = ±1 and 7 9^ the model is in 
the Ising universality class described by a c = 1/2 confor- 
mal field theory (GET). Instead, in the segment 7 = 0, 
I ft, I < 1 the underlying GET has central charge c = 1. 
To specify completely the problem we must fix boundary 
conditions (BG's). As is customary to avoid unnec- 
essary complications we will fix BG's on the fermions 
. Diagonalization brings Eq. ((l} to free Eermion form: 
H = J2k'^^kvlvk- Our choice of BG's fixes quasimo- 
menta to be quantized according to A: = (2n -I- l)7r/_L, 
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Figure 1. Typical behavior of £{t). The inset shows Gaussian 
behavior for short times, as happens for the pure case Q- 
Here L = ICQ, 13 = 6, /io,i = 1, 70 = 0.5, 71 = 0.8. 



-L/2, . . . , L/2 — 1, whereas the single-particle dis- 



persion is Afe = Y (cos k + h) -I- 7 2 sin k. 

The Loschniidt echo has been shown for the XY chain 
to be ,8] Lit) =Uk>f^fk {Alt), with 



fk (Alt) - 



l + V4-(4-l)«fcsin^(Afc^) 

1 + Ck 
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where c,- = cosh (/3A0), = sm^AOk) , A^^ = 61 - 61 
and 6k — arctan [7sin (fc) / (ft, 4- cos (fc))]. From its ex- 
plicit form we can read off a number of important points 
which we will use extensively in the following: i) the time- 
dependence is governed by L/2 frequencies A^., ii) the LE 
is a product of an extensive number of terms, and in par- 
ticular iii) the LE is a product of L/2 functions over the 
L/2 allowed values of k. The dependence on k is ana- 
lytic everywhere except for the critical points (7 = and 
< 1 or = 1 and 7 7^ 0). No singularity other than 
those expected at criticality emerges. 

Typical behavior of C{t) is depicted in Fig. [T] The 
LE quickly drops from unity at i = and then oscillates 
about its average value, with almost periodic revivals @. 

Following the spirit of Refs 0, ) we are interested in 
the distribution function of the LE seen as a random vari- 
able over infinite time equipped with the uniform mea- 
sure. The probability density of the LE can be written as 
Pc {x) := 5 {£ (t) — x), where the bar denotes the time 
average (i.e. / = liuiT^oo T^^ f (t) dt). Saying that 
the LE spends most of the time close to a certain value 
corresponds to a concentration result for Pc {x) . 

The moments of the LE can be computed us- 
ing the methods developed in Q- Here one has 
the additional complication given by the presence of 
the square-root in Eq. ([2]), which must first be ex- 
panded into an infinite scries. The result for the 
first moment is £ = nfe>o/fe' '^ith fl — ^ ~ 
(1 - ^ + [lE (6,) + 6fe/4 - 1] . Here, = 

(1 — Cj7^) sin^ (A6'fe) and E is the complete elliptic in- 
tegral of the second kind. Expanding in the small 
quench regime, that is up to second order in A0fc, one is 



able to relate the dynamical quantity C to a static quan- 
tity. Specifically, one obtains £ F (po,Pi)^, where po,i 
are Gibbs states with Hamiltonians i?o.i- This result 
extends the pure state result C = tr (p^) ~ |('0o|^i)|'^ 
which can be recovered sending /3 — >■ cxj 

The distribution function for the LE in the Ising model 
(i.e. 7 = 1) at zero temperature was considered in [3|. 
Through numerical simulations it was argued that, in 
the off-critical regime, two different behaviors were ob- 
served. The distribution of the LE was seen as similar 
to an exponential one, (Pc (x) ~ i? (x) e~^l '~ / C) or to a 
bell-shaped Gaussian-looking one. In the next section we 
will unify both of these conjectured results. 

Off- critical regime and Gaussian equilibration The 
form of the LE suggests that the LE should be thought 
of as a product of variables. Let us then consider the 
new variable Z = InC We will show that, under a 
very mild hypothesis, the variable Z satisfies the stan- 
dard central limit theorem (CLT). In particular, in the 
off-critical regime, as L — > 00, the rescaled variable 
Y = (^Z — Z) /\fL will tend in distribution to a Gaus- 
sian with zero mean and well-defined variance. To this 
aim we will show that all the cumulants of Z scale ex- 
tensively, so that for the rescaled variable Y we will get 
K„ (r) cx for n > 2 while ki {Y) = by con- 

struction. Hence only the first two cumulants of Y sur- 
vive in the L — > 00 limit, thus showing Gaussianity of 
Y . In turn, Gaussianity of Y implies that the LE is ap- 
proximately Log-Normally distributed. This explains the 
behavior observed in [tI, as a Log-Normal has regimes 
where it looks approximately exponential or Gaussian. 

In order to prove our assertion we need the (logarithm 
of the) moment generating function of Z, (A) :— 
e^^ — C^. At this point we make the reasonable as- 
sumptions that the L/2 frequencies A^ are rationally in- 
dependent (that is, linearly independent over the field 
of rational numbers). Thanks to rational independence 
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(RI) we can use the theorem of averages (see e.g. 
page 286) to compute the time-average of as a phase 
space average over an L/2-dimensional torus [3]. Our 
numerical simulations show that a possible rational de- 
pendence is very mild and it would be quite unlucky to 
produce enough correlations to invalidate the CLT. With 
RI, we obtain 

(A) = [] gk (A) , 5fe (A) = ^ [fk dd. 

k>0 ^TT Jo 

Hence (A) = exp^^^g Injifc (A). The last steps of 
the proof come from the fact that In 17^ (A) as a function 
of k is Riemann integrable, with a finite integral, provided 
we are away from critical points. Moreover, in the same 
region of parameters. In gk (A) (and so its integral over 
k) is analytic in A. Specifically, for large L, we obtain. 
In [M^ (A)] ~ LG (A), with G (A) = /J" In g^ (A) dk/ (27r) 
analytic in A. Differentiating with respect to A we obtain 
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Figure 2. Pz close to the Ising (left) and anisotropy transition 
(right). As L grows we enter the off-critical regime and Pz 
becomes Gaussian. Close to the quasi-critical regime (small 
L) the distribution becomes a broad, generally double-peaked 
function. For the anisotropy transition, one can have L for 
which the highest amplitudes are nearly equal (see text). This 
results in a collapse from two peaks to one. Parameters are 
/3 = 40 and (left) ho = 0.98, hi = 1.02, 70,1 = 1.0 and L = 50 
to 200 in steps of 30, (right) /lo.i = 0.5, 70 = 0.01, 71 = -0.01 
and L = 50 to 100 in steps of 10. Another way to enter 
the off-critical regime is to increase the temperature. Similar 
plots are obtained replacing L with the temperature T. 



that all the cumulants of Z are extensive, which com- 
pletes the proof. □ 

In particular, one has the CLT anywhere away from 
the critical points: no other source of singularity emerges 
other than those expected at criticality. 

Let us now pause for a moment and discuss how the 
CLT could be violated. One possibility is that the 
variance of Z may grow with L more than extensively, 
i.e. K2 {Z) (X L'^, with Q > 1. This would imply that 
the variance of the rescaled variable would diverge as 
L 00, thus breaking the CLT. It can be shown that 

K2 (Z) = J2k>0 ^2 (k) with K2 (k) = 1712 (fc) - [mi (fc)]^, 
and m„ (k) = ^ " [In { fk (i9))]" dt-) with n = 1, 2. By 
direct inspection of the integrals it turns out that K2 (k) is 
a bounded function in the entire parameter range. Hence 
'^2 (Z) < const. X L also at critical points. 

Quasi- critical regime and universal critical equilibra- 
tion In Ref. it was argued that for a small quench 
close to a critical point, no observable (except for trivial 
constants of motion) thermalizes. Here we will show that 
this result generalizes to the mixed case considered here. 
Moreover, as we will see, some universal features of the 
underlying critical theory show up in the long-time dis- 
tribution function. For the reasons explained above, the 
right quantity to look at is the Log of the LE. 

Since we are interested in the small quench regime, we 
expand the Log of the LE up to the first non-zero order 
in Adk- The constant terms add up to contribute to the 
average and, dropping fourth-order terms and going to 
the energy variable ujj = 2Aj.^ we arrive at 



lnC{t) ~ Z + ^ Oj cos (tujj) , 
3 

where the amplitudes are defined via a (k) 
(1 — c^^) (A0fe) and aj = a (kj). 



(3) 



Now we make the important observation that the quan- 
tity ^ is in fact a sum of i/2 independent random vari- 
ables. This can be shown assuming again RI of the fre- 
quencies ujj . Using the ergodic theorem one realizes that 
the moment generating function of ln£ is simply the 
product of L/2 generating functions. Taking the Fourier 
transform, one sees that each variable is distributed ac- 



cording to Pj {x) 



c^, with zero 



mean and variance a|/2. 

We are now in a position to understand what can 
happen at criticality and in which sense we can ex- 
pect violation of the CLT. As explained above, the to- 
tal variance, which in the small quench regime reads 
K2 (Z) — (1/2) ^^ a^, cannot grow more than exten- 
sively. But the other extreme is possible, namely the 
variances a|/2 can go to zero as L increases, and this 
can happen for most of the L/2 variables. When this is 
the case, Eq. ^ effectively represents a sum of very few 
independent variables, and the CLT regime cannot be 
reached. Here we notice that for an infinitesimal quench 
K2 is related to the fidelity susceptibility, a central ob- 
ject in the so-called fidelity approach to quantum phase 

transitions Hill- 



As we will see, close to criticality aj is a rapidly- 
decreasing function of j, so that only few amplitudes 
are appreciably different from zero. In this situation, 
a good approximation to the distribution function for 
Z is given retaining the rimax largest amplitudes aj in 
Eq. dll). Choosing rimax — 1, the distribution is the just- 
encountered Pj,-^^^ {x) with square-root singularities at 
i'^imax- With n,nax = 2 the distribution is still a very 
spread double-peaked one, with logarithmic singularities 
at Z± ||ai| — 1 02] I as shown in f?\. Using the ergodic the- 
orem it can be shown that this distribution is precisely 
the density of states (DOS) of a tight-binding model in 
two dimensions, with anisotropic couplings. In general, 
the distribution function obtained by keeping nmax am- 
plitudes is the density of states of a hypercubic rimax- 
dimensional tight binding model with anisotropic cou- 



plings aj/2 {j 



1, 



in each direction. Adding 



more and more amplitudes, eventually the CLT sets in 
and the distribution approaches a single-peaked Gaus- 
sian. Clearly, when nmax is small the distribution func- 
tion is very spread with a large variance, so thermaliza- 
tion does not take place. 

Let us now discuss the behavior of Oj close to crit- 
icality. The XY model has two different kinds of criti- 
cal regimes characterized by different underlying effective 
field theories. We now consider separately both critical 
regimes. First of all, note that increasing the tempera- 
ture simply has the effect of multiplying a (fc, T = 0) by a 
factor (1 — cosh~^ (Afc/T)) < 1. At the Ising transition 
we observe a large peak in a{k) close to k = n. The rea- 
son for the peak has to be ascribed to the single-particle 
energy vanishing as uj = v {k — tt) (where ?; = 2 I7I is a 
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velocity). The precise mechanism has been explained in 
[lO| for the pure case. At finite size the quasimomenta k 
take only discrete values. Correspondingly, most of the 
weight is absorbed by those fc's which fall in the peak. 
Other amplitudes a{kj) are considerably smaller. As a 
result, a good approximation to the distribution can be 
given by a 2D DOS as shown in Fig. [2l left panel. 

The situation at the anisotropy transition (c = 1 line) 
is very similar, with some notable difference due to the 
precise character of the c = 1 CFT. As can be easily 
seen, a (k) now has two peaks, due to the presence of 
two chiral (Majorana) Fermions corresponding to the two 
branches ofw = v\k — kp\. The doublc-pcakcd form of 
a (k) has some detectable consequence on the structure 
of the distribution function. Namely, according to differ- 
ent quantization of quasimomenta (and damping factor 
due to temperature) the allowed values of k can fall sym- 
metrically displaced among the peaks. When this is the 
case we will observe, somehow accidentally, a distribu- 
tion function given the 2D DOS with ai =02. In this 
case the two peaks of the distribution merge into a single 
one, as can be seen in Fig. [2] right panel at L = 60, 90. 

Generalization We now give an argument in sup- 
port of the validity in general of this scenario for small 
quenches. Let us restrict, for simplicity, to zero temper- 
ature. Assuming a completely generic, non-degenerate 
Hamiltonian H = En\n){n\, the LE reads £ {t) = £+ 

'^J2n>7n PnPm cos {t {En - E,n)), whcrC Pn = \{n\lpo)f 

for an initial state |V'o)- Consider now the logarithm 
of the LE and expand it in the small quench param- 
eter (that is in the perturbing potential V, which we 
assume to be extensive). Up to second order we ob- 
tain \jiC{t) = Z + 2J2n>o Pn cos {t {En- Eq)), whcrc 
for a smaU quench p„ — |(n|T^|0)|^ / {En — E^f . If we 
now assume additionally RI for the energy gaps, we re- 
turn to the previous situation with Oj = 2pj, namely 
CLT away from criticality, meaning Gaussian equilibra- 
tion. Note that the total variance is at most extensive: 
K2 {Z) = 2X;„>o P^ < 2X;„>oP« = 2x, where x is the 
fidelity susceptibility and is extensive by the extensivity 
of V and the assumption of non-criticality [l^. In the 
quasi-critical regime only a few terms of the sum domi- 
nate, thus breaking the CLT and leading to a universal, 
poorly equilibrating regime. 

Conclusions In this letter we have considered the fi- 
nite temperature generalization of the Loschmidt echo 
(LE) after a quantum quench. We have proved, un- 
der a very mild hypothesis, that away from critical 
points the LE is Log-Normally distributed, whereas for 
small quenches close to criticality the distribution ap- 
proaches that of the density of states of a D-dimensional 
anisotropic tight binding model, where D can be consid- 
ered small (e.g. D — 1,2). Although these results could 
be obtained analytically for the XY model considered 



here, we conjecture that such behavior is in fact general 
and not restricted to solvable models. 

LCV gratefully acknowledges support from Euro- 
pean project COQUIT under FET-Open grant number 
2333747, NTJ from an Oakley Fellowship, and PZ from 
NSF grants PHY-803304, DMR-0804914. 



[1] M. Rigol, V. Dunjko, and M. Olshanii, Nature 452, 854 
(2008). 

[2] J. Deutsch, Phys. Rev. A 43, 2046 (1991). 

M. Srednicki, Phys. Rev. E 50, 888 (1994). 
[3] J. von Neumann, Zeit. Fiir Phys. 57, 30 (1929), see also 

the English translation: Eur. Phys. J. H, 35, 201 (2010). 

S. Goldstein, J. L. Lebowitz, C. Mastrodonato, R. Tu- 

mulka, and N. Zanghi, Phys. Rev. E 81, 011109 (2009). 

H. Tasaki (2010), arXiv:1003.5424. 
[4] A. Uhlmann, Rep. Math. Phys 9, 273 (1976). 
[5] E. Barouch, B. M. McCoy, and M. Dresden, Phys. Rev. A 

2, 1075 (1970). 

[6] To be specific, we will use anti-periodic BC's for the 

Fermi operators, i.e. cl+i = — ci. 
[7] L. Campos Venuti and P. Zanardi, Phys. Rev. A 81, 

022113 (2010). 

[8] P. Zanardi, H. T. Quan, X. Wang, and C. P. Sun, 

Phys. Rev. A 75, 032109 (2007). 
[9] J. Happola, G. Halasz, and A. Hamma (2010), 

arXiv:1011.0380. 
[10] L. Campos Venuti and P. Zanardi, Phys. Rev. A 81, 

032113 (2010). 

[11] D. Rossini, T. Calarco, V. Giovannetti, S. Montangero, 
and R. Fazio, J. Phys. A: Math. Theor. 40, 8033 (2007). 
D. Rossini, T. Calarco, V. Giovannetti, S. Montangero, 
and R. Fazio, Phys. Rev. A 75, 032333 (2007). 

[12] V.I.Arnold, Mathematical Mettiods of Classical Mechan- 
ics (Springer- Verlag, 1989), 2"^^ ed. 

[13] Given the form of the Ak we can expect that RI holds 
for most ^,h (i.e. except for a set of zero measure) and 
for some L. From Gauss' theorem on the irreducibil- 
ity of cyclotomic polynomials (see e.g. [31 Chapter 12) 
one can derive rational independence of {cos (27rj/I/)} 
for j — 0, 1, . . . , (L — 1) /2 for L prime. Calling = 
expii9j with — 2Tij/L, Gauss' theorem asserts that 
X^jHo^ '^j Cj — with Uj £ Z implies nj = whenever L is 
prime. Taking real and imaginary parts one obtains no -f 
E-t7''^'<cos(^,) = and Y.%T'^'^ sm{-d,) = 0, 
with = n-j ± nj_|_(i_i)/2- Since the numbers nt and 
n~ are independent, the result follows. Hence one can ex- 
pect that such RI carries over to the A^, but it is possible 
that the functional dependence may lift the requirement 
that L is prime. 

[14] P. Zanardi and N. Paunkovic, Phys. Rev. E 74, 031123 
(2006). 

P. Zanardi, P. Giorda, and M. Cozzini, Phys. Rev. Lett. 
99, 100603 (2007). 
[15] L. Campos Venuti and P. Zanardi, Phys. Rev. Lett. 99, 
095701 (2007). 

[16] J.-P. Tignol, Galois' Theory of Algebraic Equations 
(World Scientific, 2001). 



